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Abstract. We propose an exactly solvable model of one-dimensional anyons with 
competing i5-function and derivative (5-function interaction potentials. The Bethe 
ansatz equations are derived in terms of the iV-particle sector for the quantum 
anyonic field model of the generalized derivative nonlinear Schrodinger equation. This 
more general anyon model exhibits richer physics than that of the recently studied 
one-dimensional model of (5-function interacting anyons. We show that the anyonic 
signature is inextricably related to the velocities of the colliding particles and the 
pairwise dynamical interaction between particles. 
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1. Introduction 

Anyons - which interpolate between bosons and fermions - may exist in two dimensions 
(2D), obeying fractional statistics. Fractional statistics have recently been observed 
in experiments on the quasiparticle excitations of a 2D electron gas in the fractional 
quantum Hall (FQH) regime |1] . The anyonic quasiparticles have profound implications 
for topological quantum states of matter [2]. For example, non-Abelian topological 
order can be studied through manipulating FQH states in a 2D electron gas. As 
a consequence, the concept of anyons has become important in topological quantum 
computation [3], H]. In one dimension (ID) anyons acquire a step-function-like phase 
when two identical particles exchange their positions. Different aspects of anyons in ID 
have been considered [5],[6l[7j. Among these an anyonic extension [8] of the ID integrable 
Bose gas with ^-function interaction [9l [10] has attracted recent attention, including the 
basic construction of the anyon model [HI [121 [13] . Many new results have been obtained 
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for this model, including the connection with Haldane [H] exclusion statistics [151 [TB] . 
correlation functions [121 [13 [HI [IS]) entanglement [20] and expanding anyonic fluids 

A key feature of ID anyons is that they retain fractional statistics in quasi- 
momentum space and the anyonic signature of the exchange interaction - the topological 
anyonic interaction and dynamical interaction are inextricably related. ID anyons show 
an intriguing sensitivity of the boundary condition of their wave function on the specific 
position of the particles [HI [121 [23| . Therefore the study of exactly solvable interacting 
anyons in ID should give insight into understanding topological effects in many-body 
physics. 

In this communication we propose a ID model of anyons with 5-function and 
derivative (5-function interaction and solve the model by means of the co-ordinate Bethe 
ansatz. The model contains an additional independent interaction parameter so that 
in principle one can study anyonic signatures through tuning the interactions. For 
this more general model we show that the anyonic signature is inextricably related to 
the velocities of the colliding particles and the pairwise dynamical interaction between 
identical particles. 



2. The model 

We consider creation and annihilation operators '^\x) and '^{x) at point x, which satisfy 
the anyonic commutation relations 

^i{x)^!\y) = e-'^'"^''^y^^\y)^ix) + 5{x - y) 



x] 



^(x)^(y) = e"'"'(^'J')^(y)$(x). (1) 

The multi-step function w{x,y) appearing in the phase factors satisfies w{x,y) = 
—w{y,x) = 1 for X > y, with w{x,x) = 0. In terms of these operators the hamiltonian 
describing anyons of atomic mass m confined in a length L is 

H= — / dx'^UxWJx) + -g dx'^UxW(x)'^(x)'^(x) 
2m Jo 2 Jo 

+ k dx[^\x)'^\x)'^{x)'^ccix) -'i!l{x)'^\x)'^{x)'^{x)] (2) 
Jo 

where we also impose periodic boundary conditions. 

The tunable parameters in the model are (i) g = fi^c/m, the dynamical interaction 
strength which drives the particles in elastic scattering through zero range contact 
interaction, where c is the coupling constant, (ii) the nonlinear dispersion coupling 
constant e, which describes a delta shift in collisions, and (iii) the statistical parameter 
K, which varies the statistics of the particles from pure Bose statistics {k = 0) to pure 
Fermi statistics (k = it). The model differs from the interacting ID Bose gas [9l [TO] 

J Further discussion of the ground state properties is given in Ref . j22j . 
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and ID interacting anyons [S] due to the presence of the nonhnear dispersion e, i.e., the 
s-wave scattering is disturbed by a delta shift in coUisions. Effectively, the ID scattering 
length is increased or decreased in individual pairwise collisions by this dispersion. It 
follows that the model contains well known models as special cases in certain limits. 
For K = 0, hamiltonian ([2]) reduces to that of the generalized nonlinear Schrodinger 
model [211 |26l |2Z|, which has been studied by means of scattering Bethe ansatz in the 
context of soliton physics. For e = 0, it describes the model of ID interacting anyons. 
For K = e = 0, it reduces to the model of ID interacting bosons, while for k = c = it 
reduces to quantum derivative nonlinear Schrodinger model [2BJ. As we shall see, the 
general model ([2]) exhibits more exotic dynamics than the two models corresponding to 
the special cases k = and e = 0. 

For convenience, we hereafter set h = 2m = 1. We also use a dimensionless coupling 
constant 7 = c/n to characterize different physical regimes of the anyon gas, where 
n = N/ L is the linear density. 



3. Bethe Ansatz 



We first present the corresponding equation of motion — i dt'^i^x, t) = [H, \E'(x, t)] via the 
time-dependent quantum fields, i.e., the nonlinear Schrodinger equation is given by 

i^t = -^^^ + ^*^^* + 4ie^^^^^. (3) 

We note that this is an anyonic version of the generalized nonlinear Schrodinger 
equation. Restricting to its A^-particle sector we can reduce the eigenvalue problem 
of the hamiltonian ([2]) to a quantum many-body problem. To this end, we define a 
Fock vacuum state \E'(x)|0) = 0. Thus we can prove that the number operator and 
momentum operator P are conserved, where 

N= [ dx^^(x)^(x) (4) 
Jo 

-L 

da;^J^^(x). (5) 







In order to properly assign the anyonic phase ) [S [la [m [12] in the domain 

Xi < X2 < . . . < xn, write the iV-particle eigenstate as 

rL 

1$)=/ dx^e-'^''^x{xi...XN)^^{xi)...^^{xN)\0) (6) 



where the wavefunction amplitude is of the form 

Xixi ...xn)= e"'^"^^><-^' "^"^'"^^ ^(''P^ ■ ■ ■ A;p7v)e'(^'""i+-+'=^^"^). (7) 

p 

Here the sum extends over all A^! permutations P. The choice of the sign of the 
anyonic phase factor in equation ([?[) can be arbitrary and may lead to different boundary 
conditions [12]. Here we prefer to a fix a boundary condition where we count the anyonic 
phase factor in equation ([7]) through the order of the particles in the scattering states 
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[TT] . By acting on the eigenstate ([5]) with hamiltonian ([2]) the eigenvalue problem for 
hamiltonian ([2]), namely 

H\^)= dx^e-''^''^HNx{xi...XN)^\xi)...^!\xN)\^) (8) 
can be reduced to solving the quantum mechanical problem 

Hmx{xi- ■ -Xn) = Ex{xi. . .xn) (9) 



where 



Hn= - Y^Q^ + g Yl "^(^^ - 

i=l * l<i<i<Af 

+ i2e ^ + (10) 

l<i<j<Af 

This anyon model turns out to be Bethe ansatz solvable with independent values 
of c, e and k. Changing the coordinates to centre of mass coordinates X = {xj + Xk)/2 
and Y = Xj — Xk leads to the eigenvalue equation (l9j) in the form 

Q2 1 92 _ 2 52 92 



dxl 2 9X2 dY^ dx% 



+ 2c6{Y) + iAe6iY)^x-Ejxi■■■x^■■.XJ...) = 0. (11) 

Integrating both sides of this equation with respect to Y from — e to +e and taking the 
limit e — gives the discontinuity condition 

i^^Xj ^Xi) Xi.Xli ■ ■ ■ yXi, Xj, . . . , X ]\f^\xj=xi+e 
(f^Xj ) Xi.Xl 5 • • • ; Xj , Xj , . . . , X^v) \xj=Xi—e 

= [2c + i2e{dx, + dxj)]x{xi, . . .,Xi,Xj, . . .,xn)\x,=x, (12) 

on the derivative of the wavefunction. This equation gives a relation between the 
coefficients A{kpi . . .kpjsf) of the form 

• • = X -kl- tc - I't] + m^^- ■ ■ ■ ■ ^^^^ 

which is the two-body scattering relation. This relation gives rise to a functional scalar 
scattering matrix which factorizes three-body scattering processes into the product 
of three two-body scattering matrices [24j. In this sense the Yang-Baxter equation 
is trivially satisfied, guaranteeing no diffraction in the scattering process [IDl El]- 
Consequently the X-particle scattering amplitude can be factorized into the products 
of two-particle amplitudes. Another way to establish the integrability of the model is to 
explicitly construct an infinite number of higher order conserved quantities, as shown 
for the ID Bose gas in Ref. [25] . 

The effective interaction strengths c' and e' appearing in the solution are given by 

^ and e' = ^— -. (14) 



cos(ft/2) cos(k/2) 
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In this model these effective coupling constants implement the transmutation between 
statistical and dynamical interactions. 

To complete the solution in terms of the Bethe ansatz we impose the periodic 
boundary condition xi^i = 0, 0:2 ... , Xn) = x{x2, ■ ■ ■ = L) on the anyonic 

wavefunction in the fundamental region < xi < X2 < ■ ■ ■ < xn L. As in the 
solution of the interacting anyon model (e = 0), this leads to the appearance of various 
phase factors [TTl [12] . The wave functions for other regions can be determined through 
application of the anyonic symmetry. The periodic boundary condition can also be 
imposed on other wavefunction arguments Xi, in each case leading to the same anyonic 
phase factor in the Bethe ansatz equations (BAE) [HI |12] . For this model, the energy 
is given by 

N 

E = Y.kl (15) 

i=l 

where the quasi-momenta kj satisfy the BAE 

ifc.L _ _MN-i) TT kj-ki + i{c' -t'{kj + h)) 

\}^k,-h-i{d-e'{k, + h)) ^ ^ 

for j = 1, . . . ,N. The total momentum is 

N 

p = "^ki = N{N -1)k/L. (17) 

i=l 

In minimizing the energy we consider k{N — 1) = u (mod 27r) in the phase factor. 

The Bethe ansatz solution provides in principle the full physics of the model. For 
the case k = 0, the BAE (fT6i) are consistent with the scattering states constructed for 
the generalized nonlinear Schrodinger model [2H [26] . It is also clear to see from the 
BAE ([IS]) that the phase shift 

^ _ kj - ke + i{c' - e'{kj + kj)) 
^'^ ~ ^ [kj -ki-i{d - e'{kj + k^)) _ 
depends essentially on the anyon parameter k, the dynamical interaction strength c and 
the nonlinear dispersion of velocity e through the effective interactions ([Hj). Among 
the competing interactions, the nonlinearity parameter e is essential in determining the 
bound states. The (5-function interaction strength c plays the role of a driving force in 
collisions. However, the anyonic parameter k determines the statistical signature of the 
particles. 



4. Preliminary analysis 

The asymptotic behaviour of the BAE f[T6l) reveals that bound states exist for c < 
and c < 2ev/L. In general, the form of the bound states appears to be extremely 
complicated. For example, for N = 2 and L{c — pe) ^ 1, we have the bound state 
{^p ± ii(c — pe)}- In this communication we focus on the repulsive regime, where c > 
and |e| < c/{2p), for which the Bethe roots are real. 
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Figure 1. Typical figure showing the quantum drift effect of the nonhnear dispersion 
e on the quasimomentum distribution function. Here k = 0, c = 0.04, L — 2 and 
N = 45. Drifts to the left and right are symmetric with respect to the sign of e. 



In the weak coupling limit 7 <^ cos(k/2), the leading term for the ground state 
energy, obtained from the BAE ( |T6|) . is 



N 



L 



2e'u/L) + -. 



(19) 



To this order, for /t = the ground state energy is independent of the nonlinear 
dispersion e and reduces to that of weakly interacting bosons. However, the 
quasimomentum distribution 

deviates from the usual semicircle law [9l [301 EI] at e = (see figured]). For e > 0, the 
high density distribution drifts to the right hand side and the front becomes steeper as 
the nonlinear dispersion e increases. For e < 0, the high density part drifts to the left 
hand side. The nature of this kind of quantum drift in quasimomentum space is due 
to the delta shift or equivalently the nonlinear dispersion of velocity, which effectively 
increases or decreases the ID scattering length in individual pairwise collisions. 

When the total momentum is non zero, i.e. p 7^ 0, the left and right drifts are 
no longer symmetric with respect to the sign of e. Here the anyonic parameter k also 
introduces quantum drifts in quasimomentum space. A particular example is illustrated 
in figure O For e > 0, the front of the distribution become very steep and the width 
becomes narrower as e increases. The energy decreases as e increases. However, for e < 0, 
the distribution becomes flatter. In this case, the energy increases as |e| increases. 
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Figure 2. Typical figure showing the quantum drift effect in the quasimomentum 
distribution function with nonlinear dispersion e and anyonic parameter n. Here 
K = 7r/100, c = 0.04, L = 2 and N — 45. The drifts is no longer symmetric with 
respect to the sign of e. 



In the strong coupling limit Lc ^ 1 , one can perform the strong coupling expansion 
with the BAE ( fT6l) . In this limit the ground state of the model becomes that of the 
Tonks-Girardeau gas. This is mainly because in the strong repulsive limit, the particles 
behave like fermions in thermalized states. As a result the effect induced by the nonlinear 
dispersion of velocity is suppressed. To obtain the explicit result for the ground state 
energy, let all quasimomenta kj shift to = kj ~ u/ L and for simplicity, choose to 
be odd. In this limit, the quasimomenta are given explicitly by 

t2 t\^2- 



{Lc' - eV)2 J 3L{Lc' - e'u) 



2N ^ 



where Uj = —{N — l)/2, . . . , {N — l)/2. This result indicates that Aj = 0. The 

corresponding ground state energy is given by 

Eo (N^ - l)7r2 / 2N 

1 + 71-, TT + T?- (22) 



A^ 3L2 V {Lc'-e'u)J L"^' 

It can be seen from the Bethe roots ( !2T|) that the flat fermion-like distribution 
becomes inclined as |e| becomes large. The distribution almost linearly increases or 
decreases depending on the direction of the current and the parameters k, and e. The 
system is strongly collective in the Tonks-Girardeau limit. 
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5. Concluding remarks 

We have presented an exactly solved model of ID anyons with general 5-function and 
derivative 5-function interaction. The Bethe ansatz solution has been obtained by means 
of the coordinate Bethe ansatz. We have seen that the anyonic signature is inextricably 
related to the nonlinear dispersion of velocity and the pairwise dynamical interaction 
between identical particles. Competing interactions among the anyonic parameter k 
and the strengths of the 5-function interaction c and the nonlinear dispersion of velocity 
e result in more subtle bound states and scattering states. Preliminary analysis of the 
Bethe Ansatz solution in the weak and strong coupling limits has revealed drifts in the 
quasimomentum distribution as a function of the nonlinear dispersion and the anyonic 
parameter. 

The quasimomentum distribution may provide a plausible way to observe anyonic 
behaviour in a general model with 5-function and derivative 5-function interaction. In 
particular, it suggests the possibility of observing anyonic behaviour via the generalized 
nonlinear Schrodinger equation, perhaps with regard to the propagation of solitons 
in nonlinear media. On the other hand, the observed quantum drifts in distribution 
may perhaps also be seen in experiments with trapped ultracold atoms in and out 
of equilibrium [321 [33]. The Bethe ansatz solution should also provide an accessible 
way to explore the dynamics in time evolution of the generalized nonlinear Schrodinger 
equation Q, just as the time-dependent dynamics of the ID boson model has recently 
been investigated [3lj. In addition, the investigation of fractional exclusion statistics 
[m Uni Uni [351 [MI [SH [M] in this generalized ID anyon gas would provide insight into 
the influence on the statistical signature of both dynamical interaction and nonlinear 
dispersion. 
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